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Abstract 

! A new scheme for non-destructive characterization of quasi-phase-matching 

o 

grating structures and temperature gradients via inverse Fourier theory us- 

0\' 

^ ! ing second-harmonic-generation experiments is proposed. We show how it is 

O 

"a 

> 



possible to retrieve the relevant information via measuring only the power in 
the generated second harmonic field, thus avoiding more complicated phase 



^ ■ measurements. The potential of the scheme is emphasized through theoretical 



and numerical investigations in the case of periodically poled lithium niobate 
bulk crystals. 
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1. Introduction 

Quasi-phase-matching or QPM is a major alternative over conventional phase-matching 
techniques in many laser applications based on frequency-conversion processes in nonlinear 
optical media (For reviews, see^* 2 -). With the maturing of QPM by periodic poling of LiNb03 
(PPLN)^^ it has become possible to produce more complicated QPM gratings simply by 
writing of the corresponding photo-lithographic mask 6 . This has led to a tremendous activity 
in engineered QPM gratings for applications in photonics. 

A proper design of the longitudinal grating structure allows, e.g., for distortion free tem- 
poral pulse compression^, broad-band phase matching 8 , multi- wavelength SHG*^, enhanced 
cascaded phase shifts, and optical diodes^ and gates^. In high power schemes, QPM soli- 
tons are known to exist 2 and longitudinal engineering can be used to tailor the soliton o 14 ! 15 
and to increase the bandwidth for their generation^- Transverse patterning can be used for 
beam-tailoringi£, broad-band SHGU 8 -, and soliton steerin g 19 i 20 . 

Furthermore, nonuniform temperature distributions in the heated material may occur, 
leading to longitudinal dependent QPM conditions. This can result in a reduced nonlinear 
conversion efficiency 2 ^ or lead to positive effects such as a reduction of the fundamental- wave 
losses in cascaded phase shift configurations 22 . 

The work presented here aims at characterizing these QPM structures in a non-destructive 
manner. The second-harmonic-generation process (SHG) has traditionally been the pre- 
ferred choice in such attempts, since it requires only one tunable laser. By working in the 
low power regime, Fourier transform theory can be applied to retrieve the grating function. 
Until now the proposed methods have required the measurement of both the phase and the 
power in the generated SH field. Such methods are difficult to realize experimentally and 
it would be much more convenient if only the power was required. This is indeed what we 
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propose here. By utilizing a mirror in the setup we show how it is possible, in principle at 
least, to obtain all relevant information via easily performed power-only measurements. 

We start by developing the general Fourier scheme with arbitrary choice of mirrors in 
section El Afterwards we apply the scheme to the concrete case of PPLN. In section 0] 
we focus on determining the grating function. This can be done at room temperature 
because we operate in the low power regime, i.e. we can neglect photo-refractive effects. 
When heated the temperature distribution along the crystal may become nonuniform and 
the determination of such temperature profiles is the subject of section |SJ Throughout the 
paper we have made an effort to ensure that the numerical simulations are in agreement 
with the parameter settings which will be encountered in the laboratory. 

2. Theoretical description of the inverse method 

We consider type-I SHG in a lossless QPM x crystal. The evolution along the propagation 
direction z of the normalized slowly varying envelopes E\ = E±(z) and E 2 = E 2 (z) at the 
fundamental pump wavelength X p (FW) and at the second harmonic wavelength X s = X p /2 
(SH), respectively, is governed by 



The wave vector mismatch, the so called phase mismatch, is introduced through j3(z) = 
47t(?t, 1 — n 2 ) I X p where n\ = n 1 (A p , T) and n 2 = n 2 (X s , T) are the wavelength and temperature 
dependent refractive indices in the crystal experienced by the FW and the SH, respectively. 
The refractive indices and hence the phase mismatch (3 can be estimated via Sellmeier fits. 




(1) 



(2) 
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For bulk PPLN the temperature-dependent Sellmeier equation^ reads 1 

where the wavelength A is measured in [A] = /im = 10 _6 m and the temperature T is 
expressed in degrees Celsius with the temperature parameter / given by 

f = (T — 24.5 °C)(T + 570.82 °C). (4) 

Variations in the second order susceptibility x^ due to the imposed grating is accounted 
for through the grating function d(z), i.e. x^ % \ z ) = d(z)xint- Xint ^ s the intrinsic strength 
of the largest second order susceptibility coefficient exploited in QPM configurations which 
for lithium niobate is xlnt = ^33 ~ SOpm/V. Furthermore, in experiments the crystals are 
often heated to eliminate photo-refractive effects. The heating is in general nonuniform and 
therefore the phase mismatch becomes ^-dependent through T = T(z). 

The real and measurable powers P\ and P 2 at the FW and at the SH, respectively, are 
via the normalization given by 

P 1 = V A\E 1 \ 2 , P 2 = 2 V A\E 2 \ 2 , (5) 

where 77(A P ,T ) = 47r 2^ d a ■ Vo is the impedance of free space and A is the cross area of the 
Gaussian pump beam. We emphasize that for normal temperature variations, i.e. ±10°C, 
the variations in the refractive indices are small. Hence i] can be assumed to depend only 
on A p and a reference temperature T which we take to be the temperature at the beginning 
of the crystal. We also note that because of the way we have chosen to normalize the 
system we measure r] in Watts, i.e. [rf\ = W, and Ej in reciprocal meters, i.e. [Ej] = m _1 . 



1 For extraordinary polarized electric fields leading to the use of (£33: ai = 5.35583, a-i — 0.100473, as = 
0.20692, a 4 = 100, a 5 = 11.34927, a 6 = 1.5334- 2 , b x = 4.629- 7 , b 2 = 3.862- 8 , 63 = -0.89 x 10~ 8 , and 
64 = 2.657 x 10~ 5 . 



This normalization allows us to operate with the relevant physical parameters, i.e. the 
wavelengths, the grating period, the phase mismatch, and the temperature, while preserving 
a simple structure in the governing equations. 

As described in the introduction, we aim at solving the inverse problem, i.e. to determine 
the grating function d(z), knowing only the input powers at both wavelengths and the 
output power of the SH. Solving Eq. ([TJ2J) for d(z) in principle requires information on 
both the input and output amplitude at both wavelengths and the relative phases. The 
problem is significantly reduced by assuming that the FW is undepleted in which case the 
solution can be found by applying the inverse Fourier transformation on Eq. (J2J). To simplify 
the argument we keep the temperature uniform for the moment and launch only FW. The 
normalized SH intensity at the output of the crystal is then given by 

\E 2 ((3,L)\ 2 = \EZ(0)\ 2 F{d(z)}F*{d(z)}, (6) 

where d(z) = u(z)d(z) is the grating function modified with the window function u{z) which 
is 1 if < z < L and otherwise, with L being the length of the crystal. The complex 
Fourier transform has been defined as = /f^ f(z) exp(— if3z)dz and T*{f{z)} 

denotes its conjugated. In general we do not know any of the symmetry properties of the 
grating function, except that it is real, and hence we cannot hope to determine the grating 
function from Eq.© knowing only the power spectrum of the SH alone; additional phase 
information is required. However, we observe that if we could force the grating function to 
be even then by virtue of symmetry properties the Fourier transform would become real and 
in principle Eq.fjUJ) could then be solved without the phase information. 

One way in which we can force the grating function to be even, is to put a mirror directly 
at the output surface of the QPM crystal. Mathematically speaking we simply prolong the 
original grating function with its mirror image and let this new grating function of length 
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Fig. 1. Illustration of how to generate an even grating function with the help of a mirror. 

2L enter as d(z) = d(—z) in Eq. ({%]). Naturally any temperature distribution would likewise 
become an even function in z. In Fig. ^we have sketched the situation with the emphasis 
on the point that d(z) is now both even and real; the prerequisite for solving the problem 
via inverse Fourier transformation. 

Before Eq. (J2J) is integrated in the general case with a grating and nonuniform tempera- 
ture profile, we need to elaborate a bit more on the mirror. In general the waves are affected 
in two ways by the mirror: they experience a phase shift and they are attenuated. Though 
we shall later focus on metallic mirrors where both wavelengths experience a 7r phase shift, 
we here set up more general mirror conditions since other mirrors such as dielectric ones 
might be of interest. At the mirror we have 

Ej(z = 0+) = Tj expiim^Ejiz = 0~), j = 1, 2, (7) 

where Tj and rrij are the amplitudes and the phase shifts, respectively, of the reflection 
coefficients. The reflection coefficients are in general wavelength dependent. 

We can now perform the integration of Eq.(j2J in the general case of a grating function 
d(z) and a nonuniform temperature profile. The result is 

E 2 (L) = r 2 e im2 E 2 (-L) 
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+i{rle i2mi + r 2 e im *\E x {-Lf [ d(z) [cos B cos((3 z) - smBsm((3 z)} dz 

Jo 

- [r 2 ie i2mi - ^e^E^-Lf [ d(z) [cos B sin (f3 z) + sin B cos(f3 z)} dz. (8) 

Jo 

For simplicity we have set B(z) = f Q z 5(3(z') dz' where 8(3{z) = 4n[5ni(t(z)) — 5n2(t(z))]/X p 
is the ^-dependent part of the phase mismatch /3(A P , T) = (3q(X p ,Tq) + 8/3(t(z)) with the 
reference temperature T at the beginning of the crystal and total temperature T(z) = 
T + t(z) at the coordinate z. Henceforth we shall refer to t(z) as the temperature profile. 
Srij(t(z)) is the ^-dependent part of the total refractive index rij(T , X p , t(z)) = noj(T , X p ) + 
Srij(t(z)) determined from Eq. 

Eq.fjHJ) is the starting point for all further analysis in this paper and in the next section 
we show how it applies to the case of a metallic mirror. Here, however, we feel that we must 
supply a few general comments on the structure of the integrals in Eq.fjHj) which eventually 
will be solved by applying Fourier theory. The way we have set up Eq.(jHJ) indicates that 
we have chosen the ^-independent part of the phase mismatch, f3o, and the coordinate z 
to be our Fourier variables. In an experimental situation we have the option of changing 
Po through either of the two independent variables, i.e. through the pump wavelength X p 
or through the reference temperature T . It is clear that application of Fourier's integral 
theorem to any of the integrals in Eq. (JHJ) requires that the function B must be independent 
of j3o or more specifically: B must be independent of the independent variable through 
which we change j3 . Since B depends explicitly on A p through the way we have defined 
6(3 above, we can immediately rule out to determine temperature profiles by changing /? 
through the wavelength. Hence temperature profiles must be determined by changing fio 
through temperature. In theory this requires that 5rij(t(z)) be independent of T and in 
section El we show that for small variations in temperature this can indeed be assumed to be 
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so. The same assumption has been applied on waveguides in PPLN^ and thus the theory 
presented here is applicable for other materials and experimental setups. 

Once the mirror type has been decided on Eq.fjSJ) is readily solved even if it looks com- 
plicated. We notice that by proper engineering of the mirror we can make either the sum or 
the difference between the two reflection amplitude coefficients in Eq.fjHJ) vanish. It is this 
fact that allows us to simplify Eq.fjHj). Put differently we can state that because of the spa- 
tial symmetries of the mirror-expanded grating function, which forces the involved Fourier 
integrals to be real, it is essentially no longer necessary to obtain the phase information that 
we need in order to solve Eq. (JHJ). 

3. Modeling with a metallic mirror in bulk lithium niobate 

We now focus on exploiting metallic mirrors. These mirrors are relatively cheap and can 
be produced with reflection amplitude coefficients very close to 1. Though the coefficients 
are wavelength dependent, it is reasonable to assume that T\ ~ r 2 and that they vary little 
around each of the wavelengths. With both waves experiencing a mj = tt phase shift, the 
sum between the reflection amplitude coefficients in Eq. (jSJ) vanishes and the unseeded 
normalized SH output intensity as found from Eq. (jHJ) becomes 

\E 2 (P ,L)\ 2 ^ (r 2 1+ r 2 ) 2 E 1 (-L) i f 2 (p ), (9) 

where f(f3o) = Jo d(z) [sin B cos(f3oz) + cosB sin(/3oz)} dz. From Fourier's integral theorem 
we then know that 

d(z)^B= Yl^'fjf] rf M sin5= 2 2 ^<*f» , (10) 
[r{ + r 2 )Ei{-Ly (rf + r 2 )E 1 (-L) 2 

where Tj^ and are the imaginary and real part, respectively, of the inverse Fourier 
transform. Again d(z) is the grating function modified with the window function u(z) 



according to d(z) 



u(z)d(z) where u{z) = lifO<2:<L and otherwise. With no 



temperature profile present, B = 0, and f(Po) becomes an odd function in /3 . Hence we 
can apply the inverse Fourier sine transform, •7 7 S ~ 1 {/(A))} = f Jo°° f(Po) ^(Po^d-Po, an d 
determine the grating function through 



P 2j o and P\,i are the measured SH output power and the FW input power, respectively, as 
given by Eq. (JSJ). We notice that both rj and P 2y o ar e functions of (3q and that the inverse 
sine transform must be applied to the product of the two functions. 

In principle we can now determine any grating function in a uniform temperature distri- 
bution simply by measuring the SH output power as a function of the pump wavelength; 
additional phase information is no longer required. The attentive reader will have noticed 
that the arguments of the inverse Fourier transforms in Eq. (jlOj) should have been functions 
of the absolute value of the SH field \E 2 (L)\ instead of just the real value E 2 (L). Of course 
l-E^C^)! is what we measure but we need the sign information to correctly reconstruct the 
grating function and temperature profile. Luckily the sign information is easily retrieved 
by exploiting that the derivative of E 2 (L) with respect to (3q must be a smooth function 
in /3 whereas the derivative of \E 2 (L)\ is not. Whenever the sign of E 2 (L) changes we get 
non-analytical points in the derivative of \E 2 (L)\ with respect to (3q. These points can be 
located numerically and hence E 2 (L) can be determined. 

4. Simulation results: Determining the grating function 

In the following we present simulation results validating Eq. (jllj) and we discuss the de- 
termination of the grating function in uniform temperature profiles. As mentioned above 
we are now, because of the mirror, able to reconstruct any grating function simply by mea- 




(11) 
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suring the SH output power P 2 ,o as a function of pump wavelength. In this section we 
shall first verify Eq. with simulations on a perfectly periodic QPM crystal, i.e. we 

show that the contribution from the first integral in Eq. (jSJ) is negligible even if we have 
considerable mirror losses in the setup. Secondly we shall focus on the case of a perfectly 
periodic crystal with a duty-cycle different from D = 0.5. This example clearly illustrates 
the differences between the same experiments made with and without the mirror and hence 
emphasizes the strength of the presented scheme. The last part of this section is dedicated 
to a discussion of the limitations of the presented scheme. We remark that these limitations 
are not a consequence of not making phase measurements but rather an inherent problem of 
applying Fourier theory to the particular case of determining QPM gratings in periodically 
poled materials. 

We leave the determination of temperature gradients to the next section and set 
T(z) = 24.5 in Eq. Q in the following. With a pump laser tunable in the inter- 
val X p G [1.5/im 1.6/im] this yields phase mismatches in the approximate range flo G 
[—3.6 x 10 5 m _1 —3.0 x 10 5 m -1 ]. In a traditional domain inverted QPM crystal this cor- 
responds to a domain length of ~ 10/mi or around 10 000 domains in a 1cm long crystal. 
Possibly the crystals we want to characterize are several centimeters long. With a beam 
diameter of 350/im the cross area of the Gaussian pump beam is A = 9.62 • 10 _8 m 2 and the 
wavefront can be assumed planar for approximately 10cm. With a pump power of lOmW 
this yields a FW intensity of I\ ~ 100 000W/m 2 . To establish whether this intensity is in 
the low-depletion regime we use the analytical solution at perfect phase matching for SHG 
with first order QPM2£. For propagation distances less than 9cm we get that the FW is 
unaltered to within 1 x 10~ 4 of its initial value. The numerical results confirm that this is 
indeed in the low depletion regime. Regarding the detection of the SH output we shall here 
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assume that we can measure O.lnW and numerically put all values below that to zero. In 
the following we refer to this value as the cut-off on the SH power measurement. 
A. Determination of domain length in perfectly periodic crystal 

To illustrate the method and to verify it in the presence of mirror losses we here consider 
the case of determining the domain length in a perfectly periodic crystal with a duty-cycle 
D = 0.5, i.e. the domain length is the same for the d = — 1 and d = 1 domains. The grating 
function can be expanded according to 

i(z) \ 1, 0<*<A = 4 ^ sintnWA) {u) 
1-1, A<z<2A **-W.- " 

Choosing first order QPM, i.e. n = 1 in Eq. (|12p. we can rewrite Eq. (jUJ) and express the 
output SH power as 

„ 8(rl + r 2 ) 2 L 2 



-P^sinc 2 



(13) 



7T 2 Ar/ 

where sinc(x) = sinx/x. In Fig. El we show the tuning curve found from a numerical 
experiment on a 1cm long crystal with a domain length of 1 x 10~ 5 m. In the simulation 
we have used reflection amplitude coefficients of r\ = 0.95 and r 2 = 0.75 for the FW and 
the SH, respectively. The reflection amplitude coefficients have been chosen so that the 
difference between them is sufficiently big to encompass any differences we imagine could be 
encountered in the laboratory. Finally we assume that the laser shares the above discussed 
characteristic with a FW power of Pi = lOmW and scan through the pump wavelength 
interval with a step-length of AA P = 2 x 10~ 10 m which yields a total of 500 measurements. 

It is well known 2 ^ that we would also get a sinc-shaped tuning curve from the correspond- 
ing no-mirror case and that the phase information likewise is not needed to reconstruct the 
actual grating function. The numerical experiment in Fig. Elis non-trivial because it shows 
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Fig. 2. SHG tuning curve for perfectly periodic QPM crystal of length L = 1cm and with 
domain length A = 10/im. The reflection amplitude coefficients of n = 0.95 and r2 = 0.75. 
The curve is in fact discrete which is indicated with dots in the enlarged part of the curve 
shown in the inset. 
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Fig. 3. Reconstructed grating function from the SHG tuning curve in Fig. |2l The entire 
grating function is shown to the left. Two enlarged parts of the grating is shown to the 
right: The middle part in the top and the part around the end of the crystal in the bottom. 

that we can apply the approximation given in Eq. (0), i.e. that the scheme works even for 
considerable mirror losses. In fact, we have put the theoretical curve as found from Eq. (|T3j) 
on top of the curve in Fig. El and the experimental points fall on the curve to within the 
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precision inherent to the Fourier transform and to within the precision owing to the cut-off 
on the SH power measurement. 

To recover the grating function we need to generate the full tuning curve spanning the 
entire fio-axis, thus covering all the QPM peaks owing to the different QPM orders. Numer- 
ically speaking, however, we only need to take into account the first few peaks and here we 
have chosen to work with the 7 first higher order peaks, making the computational effort 
tolerable. The expansion of the tuning curve is fairly trivial in the case of the perfectly 
periodic duty-cycle D = 0.5 crystal. First we locate the maximum of the tuning curve on 
the /3 -axis. Once we know the location j3o tmax of the maximum, then Eq. (JT3|) tells us that 
the n'th order peak is located at (3q = nf3o,max- We remark that the /3o-axis and the A p -axis 
in Fig. 121 are connected through the Sellmeier Eq. Q and we notice how the absolute value 
of Pq here is a decreasing function in X p . Hence on the A p -axis the higher order QPM peaks 
are located closer to A p = whereas they are located at numerically higher values of (3q 
than the maximum value on the /3o-axis. The series of experiments is performed keeping the 
step-length AA P constant. Since the Sellmeier equation is nonlinear the step- length Af3$ on 
the /?o-axis will not be constant and numerical routines must be supplied to fix this problem. 
We can estimate the resolution in the z-domain, i.e. the inverse Fourier domain, through 
Az = 2n/(Ap x /3o }( nsc), where Po^disc are the total number of experiments. By using the 
shortest step-length we find on the original /3o-axis we get an estimate on the coarsest reso- 
lution we can expect yielding Az ~ 1.6/xm which is good enough to resolve domain lengths 
on the order of 10/im. 

We are now ready to apply Eq. (fTCj|) and perform the inverse Fourier transform. The 
result is shown in Fig. |3J The figure shows both the complete recovered function and two 
enhancements: one of the middle part and one of the part around the end of the crystal. 
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The figure showing the complete recovered function is entirely black owing to the fact that 
the value of d(z) changes 10 000 times in the interval z G [0 1cm]. From the enhanced part 
of the grating function from the middle of the crystal we observe that the domain length 
is correctly retrieved. The small oscillations around d(z) = ±1 are due to the inherent 
difficulties in making numerical Fourier transforms and can only be reduced by enhancing 
the resolution. On the other hand the oscillations are not critical in the case of PPLN since 
we now that d(z) = ±1 and cannot take any values in between. Somewhat more critical is 
the fact that we according to the close up of the part around the end of the crystal do not 
get a nice sharp cut at z — 1cm. Instead we get a slow decrease in the value of d(z) which 
eventually becomes zero as can be seen from the figure showing the entire grating function. 
We have verified that the problems in correctly retrieving the grating function around the 
end of the crystal are due to the cut-off on the SH power measurement. With a lower cut-off 
the grating function goes to zero within very few oscillations. 

In conclusion we can say that with no cut-off on the SH power measurement and with 
indefinite resolution the grating would off course have been perfectly retrieved. On the 
other hand this is what we would expect from experiments like this and as such it is not a 
consequence of the scheme we propose, i.e. experiments made with no mirror by measuring 
the phase instead would be subject to the same difficulties. Thus we have verified that 
the scheme is applicable even with mirror losses, or more precisely: with a considerable 
difference between the loss experienced in the FW and that experienced in the SH. 
B. duty-cycle errors 

We now turn towards an example which demonstrates how our scheme immediately renders 
information which would otherwise require phase measurements. In the previous part we 
investigated the perfectly periodic QPM crystal with a duty-cycle D = 0.5. Here we shall 

14 



still consider only perfectly periodic crystals but with arbitrary duty cycle, D £ [0.5 1]. 
Such a grating function can be expanded in the Fourier series 

1, < z < 2DA 



d(z) 



-1, 2DA < z < 2A 

/l ^\ 4 ^ sin(n7r/2) sm(nDn) 
- 2 (n-D) + - £ —^-^ -sm 



2 / ^n=tt,— U 

4 ^ sin((n + l)vr/2) sin(n£>7r) . 



n—z + nir D 

A V2 



n—z + mr D 

A V2 



(14) 



The problem can still be solved analytically and for the SH output at the first order QPM 
peak, n = 1 in Eq. (JHJ), Eq. (0) reduces to 



(rl + r 2 ) 2 L 2 2 . 2 [cos [q - vrD] - cos (ttD) 

—r^ sm {L>7i) 



(15) 



7T 2 A// 

where for simplicity we have set q = (j^ + /3 J ^- We observe that for D = 0.5 Eq. (|T5|) 
reduces to Eq. l[I3J). For D 7^ 0.5 we see that the solution is scaled by a factor sin 2 (D7r). 
The solution is also no longer sine shaped as compared to Eq. (fT3|) and this is an important 
observation. In a normal setup, i.e. without utilizing a mirror, the solution would indeed 
still have been sine shaped^ making it difficult to distinguish from the duty-cycle D = 0.5 
solution. Of course it would still be scaled by the factor sin 2 (D7r) but since we would like to 
determine both the duty cycle, D, and the domain length, A, additional phase information 
would still be required. In our scheme, however, this is not so. 

In Fig. 0]we have shown the outcome of an experiment where the duty-cycle D 7^ 0.5. 
All other parameters are the same as the ones used to produce Fig. El If we compare Fig. 
121 and Fig. 0] there is an obvious difference in the shape around the central peaks. Suppose 
we did not know anything about the crystal in front of us, then because of the asymmetric 
tuning curve in Fig. 0]we can immediately conclude that either the crystal is not perfectly 
periodic or the duty-cycle D 7^ 0.5. Since we do not know the tuning curve throughout 
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the entire /3o-axis, we have to make assumptions in order to further characterize the grating 
function responsible for the tuning curve. Duty-cycle errors are not uncommon and hence it 
would be a reasonable starting point to assume that we are dealing with exactly that, then 
perform the analysis, and then verify the nature of the error by simply holding it up against 
the theoretical solution. In order to determine the duty-cycle and the domain length we 
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Fig. 4. SHG tuning curve for periodic QPM crystal with duty-cycle D = 0.7. The other 
parameters are the same as for Fig. El 

need to measure two things in Fig. 0] the maximum value, P2,max, of the generated SH and 
its location, /3o, ma z, on the /3o-axis. The maximum value Pi,max is independent of where on 
the /?o-axis it is located and hence we can determine the duty-cycle as a function of p2,max- 

Knowing D we can also find q max which is determined through the transcendental equation 

sin (f ) sin (f - nD) 
Q= sin( q -nD) • (16) 

Since we have measured f3 0ymax and q max = + Aj.max) L we can now easily retrieve the 

domain length A. In Fig. 03 we have sketched how we determine the duty cycle and domain 

length with the letters indicating the order discussed above. The point at 'A' indicates 

P^max which we have measured on Fig. |U to be P2, ma x ~ 78nW =>■ p ~ 0.6 where p = 

P2,max/P2,max(P > = 0.5), i.e. p is the maximum output power relative to the maximum of 
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Fig. El Following the line to point 'B' we determine the duty-cycle to be D m 0.70. We then 
jump to point 'C on the curve showing the domain length as a function of the duty-cycle. 
Following the line to point 'D' finally gives us the domain length A « 10.000/xm. 

The method outlined above does not involve Fourier transformation since we have already 
made assumptions as to the nature of the grating function. To verify if we are indeed dealing 
with duty-cycle errors and that we have found the right values for the duty-cycle and the 
domain length we only have to put the analytical solution (fTK|) on top of the curve in Fig. 
0] and see if it fits. We emphasize that without using the mirror setup, it is necessary also 
to measure the phase of the SH in order to apply the method illustrated by Fig. El 

110.002 



Q. 0.5 




9.992 



Fig. 5. Determination of duty-cycle D and domain length A in the presence of duty-cycle 
errors. The full curve shows the maximum relative SH output power p = P2,max/P2,max(D = 
0.5) as a function of D. The dashed line shows A, also as function of D, determined via Eq. 
ITfij) . See text for explanation of the letters A, B, C, and D. 



C. Discussion on the limitations of the inverse method 

We shall keep this discussion on the limitations of the inverse method short since these are 
a consequence of the nature of the Fourier transform rather than due to the here proposed 
scheme with the mirror setup. 
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Essentially the biggest hurdle for the inverse method is the limited /3o-bandwidth owing 
to the limited wavelength range for the used lasers. This is of course connected to the fact 
that using PPLN moves the central peak on the tuning curve far away from /3 = which 
means that the /3o-bandwidth in practice must be huge to encompass all the information 
of the SH tuning curve. Even if we had an unlimited wavelength range it is doubtful if we 
could go all the way to /3q = 0. This is evident from for instance the Sellmeier Eq. (J3J) we 
here apply for bulk PPLN. 

With a limited /3o-bandwidth far away from zero we can in theory only hope to char- 
acterize a limited category of grating functions, namely those which are periodic and for 
which the corresponding Fourier frequencies falls within the /3o-bandwidth. Aperiodic func- 
tions in general have dense Fourier spectra and hence we would not be able to get all the 
necessary information to retrieve the grating function. In a traditional setup with no mirror 
grating functions with stochastic boundary errors, with missing domain reversals or with 
single domains of different length like the optical diod©^ are all known^ to share the familiar 
sinc-shape and as such they are indistinguishable from the perfectly periodic crystal. Only 
when the relative sizes of the errors become sufficiently large compared to the domain length, 
they change the shape of the sine significantly. The situation does not change even if we use 
the mirror setup. The big problem is that the information from these aperiodic errors is not 
folded around the QPM peaks but rather around (3q = 0, far away from any realizable (3q. 

Once again we emphasize that the above limitations will always be present, mirror or 
not. However, as exemplified above for duty-cycle errors, the mirror setup allows for power- 
measurement-only characterization of periodic functions which are bandwidth limited to 
within the realizable /?o-interval. We have verified that the scheme can also be applied to 
the case of a linear increase in the domain length. 
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5. Simulation results: Measuring temperature profiles 

We now turn towards experimental conditions where a temperature profile, t(z), is present. 
In Fig. (BJwe have sketched such a situation emphasizing how the temperature profile, like 
the grating function, becomes an even function in z when utilizing the mirror setup. As 
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Fig. 6. QPM crystal in mirror setup illustrating the even mirror-expanded temperature 
profile, t(z). 

discussed in section |2] we must change /3 through varying the temperature rather than the 
wavelength in order to determine temperature profiles. Deriving the temperature profile 
from Eq. is not trivial because the refractive index through the Sellmeier Eq. (JHJ) 

is a nonlinear function in temperature. The implications of the Sellmeier equation being 
nonlinear are two-fold. First of all the function B = J 2 5j3(z')dz' must be independent of the 
reference temperature To. In theory this must hold true everywhere on the /?o-axis, which of 
course is impossible to assure. In practice however, it suffices that it holds true only within 
the narrow intervals on the /?o-axis where the non-trivial part of the tuning curve most often 
is located. In the following we verify through simulations that this is indeed so. The second 
implication of the Sellmeier equation being nonlinear concerns the last step of retrieving 
the temperature profile since this involves determining the temperature at a propagation 
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distance z knowing the refractive index at that coordinate. Assuming that the variations in 
temperature owing to the temperature profile are small, we can retrieve the temperature by 
Taylor expanding the Sellmeier equation. 

When B(z) ^ the tuning curve, essentially /(/3o) in Eq- Q, is no longer an odd function 
in (3q. Hence the expansion of the tuning curve to the negative /?o-axis is more complicated 
than before. However, here we shall again consider the case of the perfectly periodic square 
grating for which the expansion is reasonable straight forward. The qualitative shape of the 
tuning curve is determined by the function f(flo) which becomes 
2 r L 



dz, n = 1,3,. ..(17) 



The integral of the first cosine gives us the structure of the /?o > part of the tuning curve 
and the second integral the /3 < part. Since B(z) is assumed independent of /3 , this term 
distorts the QPM peaks on the fio < 0-axis in the same way as they do on the (3$ > 0-axis 
and consequently f(nn/A + X) = —/(—nir/A + X) where X is some displacement from the 
n'th peak. 

From Eq. (|l(Jj) we get that B(z) can be determined by 

BW = A r ctan(^Sl (18) 

where ^Jr]P 2) o according to the discussion above must be expanded over all the /3 -axis. 
Knowing B(z), it is trivial to retrieve the temperature profile. In Fig. [7|we show the tuning 
curve resulting from a numerical simulation on a 5cm long perfectly periodic crystal with a 
temperature profile 

t z , z G [0 0.5cm] 
*(*) = { 0.5cm • t , ze [0.5cm 4cm] (19) 
—t z/2 , zG [4cm 5cm]. 
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Besides the fact that we have used a longer crystal, the simulation has been done with 
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Fig. 7. SHG tuning curve for periodic QPM crystal with a temperature profile varying 
according to Eq. (fTHj) . The parameters are: L = 5cm, to = 1000°C/m, D = 0.5, r\ = 0.95, 
T2 = 0.75, A = 10/im, and X p = 1.6yum. The dashed vertical line indicates the location of 
the maximum with B(z) = 0. 

exactly the same numbers as used in the previous sections including the mirror losses. The 
only difference is that we have now scanned through the reference temperature To, with 
step-length AT = 0.04°C, instead of through the pump wavelength which we have fixed 
at X p = 1.6/im. With no temperature profile, the sample is phase matched at \Pq\ = n/A 
which is reached at T ~ 115°C The presence of a temperature profile changes the phase 
matching condition and we observe that the tuning curve on Fig. [7|is no longer sinc-shaped 
and that the main peak is shifted towards a higher absolute value of (3$. As we would expect 
from Eq. ()13j) . we also observe that the maximum measured powers are around 25 times 
higher than the maximum powers measured in the experiments on the 1cm long crystal in 
the last section. 

In Fig. |H1 we have plotted the retrieved temperature profile found from Fig. via Eq. 
f!18|) together with the retrieved profile from a numerical experiment with to = 500°C/m. 



21 



We see that we have excellent agreement between the actual profile and the retrieved profile 
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Fig. 8. Retrieved temperature profiles for two different t : (A) t = 1000°C/m and (B) 
to = 500°C/m. The full curves are the retrieved profiles and the dashed curves are the 
actual profiles used in the numerical experiments. 

and that the agreement seems to be better for the to = 1000°C/m simulations. This is 
because we have used a larger reference temperature interval, T G [30°C 205°C], for the 
t = 1000°C/m curve than for the t = 500°C/m curve for which T G [100°C 125°C]. We 
have used the larger interval to illustrate that B(z) can in fact be assumed independent of 
/3q over a considerable interval. The large oscillatory behavior just before z = 5cm is a result 
of the cut-off on the SH power measurement, i.e. the oscillations are not present if we lower 
the cut-off. t(z) will of course not remain independent of the reference temperature over an 
interval this large. To illustrate that we still gain information with a narrower interval we 
show the to = 500°C/m simulation. For longer crystals the tuning curves will in general be 
even narrower. 

We remark that temperature profiles like the ones depicted in Fig. |H] have been de- 
termined before from experimental measurements on uniform lithium niobate^. However, 
those experiments relied on the oven producing an even temperature profile, i.e. the ab- 
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solute values of the linear increase and decrease at the beginning and at the end of the 
oven, respectively, were the same. The profiles from Fig. |H] are odd and hence cannot be 
determined through Fourier analysis without either aquiring phase information or, as we 
have done here, by utilizing the mirror setup. 

6. Conclusion 

In summary we have investigated how the mirror-setup paves the way for characterization of 
grating functions and temperature profiles via second-harmonic-power measurements only, 
i.e. without additional phase information. The grating function and temperature profile 
both become even functions in the propagation coordinate because of the mirror and we 
derived the Fourier scheme, including losses and phase shifts due to the mirror, to solve 
the inverse problem, i.e. to find the grating function and temperature profile knowing only 
the second harmonic tuning curve as a function of either wavelength or temperature. We 
verified the scheme through numerical simulations on bulk PPLN and showed how to retrieve 
information which is bandwidth limited to within the realizable phase-mismatch interval. 
In particular we investigated the case of the perfectly periodic duty-cycle waveguide for 
which we can determine both the domain length and the duty-cycle by simply looking at 
the tuning curve resulting from the mirror setup. Concerning temperature profiles we found 
that if we know the grating function then in theory we can determine any temperature 
profile. In particular we investigated the often encountered step-profile and we saw how this 
was beautifully retrieved. 
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